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Abstract 



We present a technique for implementing in a fast way, and without any approxima- 
tions, higher-order calculations of partonic cross sections into global analyses of parton 
distribution functions. The approach, which is set up in Mellin-moment space, is partic- 
ularly suited for analyses of future data from polarized proton-proton collisions, but not 
limited to this case. The usefulness and practicability of this method is demonstrated for 
the semi-inclusive production of hadrons in deep-inelastic scattering and the transverse 
momentum distribution of "prompt" photons in pp collisions, and a case study for a future 
global analysis of polarized parton densities is presented. 



1 Introduction 

High-energy spin physics has been going through a period of great popularity and rapid de- 
velopments ever since the measurement of the proton's spin-dependent deep-inelastic structure 
function g^ by the EMC more than a decade ago. As a result of combined experimental and 
theoretical efforts, we have gained some fairly precise information concerning, for example, the 
total quark spin contribution to the nucleon spin. Yet, many other interesting and important 
questions, most of which came up in the wake of the EMC measurement, remain unanswered 
so far, the most prominent "unknown" being the nucleon's spin-dependent gluon density, Ag. 
Also, polarized inclusive deep-inelastic scattering (DIS) data do not provide enough informa- 
tion for a complete separation of the distributions for the different quark and anti-quark flavors 
u, u, d, d, s, and s. Here in particular a possible flavor asymmetry in the nucleon's light sea, 
Au — Ad ^ 0, has attracted quite some interest, and several models have been proposed recently 
0, |3], |J. Current and future dedicated spin experiments are expected to vastly broaden our 
understanding of the nucleon spin structure by studying reactions that give further access to 
its spin-dependent parton distributions, among them Ag and Au, Ad. In addition to lepton- 
nucleon scattering, there will also be for the first time information coming from very inelastic 
polarized pp collisions at the BNL Relativistic Heavy-Ion Collider RHIC || . 

Having available at some point in the near future spin data on various different reactions, one 
needs to tackle the question of how to determine the polarized parton densities from the data. 
Of course, this problem is not at all new: in the unpolarized case, several groups perform such 
"global analyses" of the plethora of data available there || |7j . The strategy is in principle clear: 
an ansatz for the parton distributions at some initial scale fio, given in terms of appropriate 
functional forms with a set of free parameters, is evolved to a scale /j,f relevant for a certain data 
point for a certain cross section. Then the parton densities at scale lip are used to compute the 
theoretical prediction for the cross section, and a \ 2 value is assigned that represents the quality 
of the comparison to the experimental point. This is done for all data points to be included 
in the analysis, and subsequently the parameters in the ansatz for the parton distribution 
functions are varied, until eventually a minimum in \ 2 is reached. 

In practice, this approach is not fully viable if the partonic scattering is treated beyond the 
lowest order of perturbation theory. The numerical evaluation of the hadronic cross section at 
higher orders is usually a rather time-consuming procedure as it often requires several tedious 
numerical integrations, not only for the convolutions with the parton densities, but also for 
the phase space integrations in the partonic cross section. The fitting procedure outlined 
above, on the other hand, usually requires thousands of computations of the cross section for 
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any given data point, and so the computing time required for a fit easily becomes excessive 
even on modern workstations. We note that for practically all reactions of interest in the 
unpolarized and polarized cases the first-order QCD corrections to the respective partonic 
cross sections are known by now. They are generally indispensable in order to arrive at a 
firmer theoretical prediction for hadronic cross sections; for instance the dependence on the 
unphysical factorization and renormalization scales is reduced when going to higher orders in 
the perturbative expansion. Only then can one reliably extract information on the parton 
distribution functions. 

In the unpolarized case, a way to get around this problem is based on the fact that the 
parton densities are already known here rather accurately || [/]]. Their gross features are 
basically determined by the wealth of very precise DIS data which cover a wide kinematical 
range in the momentum fraction x and the scale /j,p ~ Q. As a consequence, the theory answer 
for a certain cross section is expected to change in a very predictable way when going from the 
lowest-order Born level to the first-order approximation. It is then possible to pre-calculate a 
set of correction factors Ki (i running over the data points), and to simply multiply them in 
each step of the fitting procedure to the lowest-order approximation for the cross section, the 
latter being usually much faster to evaluate than that involving higher order terms. The Ki 
usually hardly change at all from one set of parton distributions to another, and in any case 
one may update them if necessary at certain stages of the fitting procedure. 

It should be noted, however, that this way of treating next-to-leading order (NLO) cor- 
rections in a fit is not necessarily adequate in all cases of interest. In particular if one is 
interested in extracting information about the gluon density at large values of x, where it is 
only rather poorly constrained at the moment, the correction factors Ki cannot be reliably 
pre-calculated, and they may vary considerably during the fitting procedure. It is therefore 
desirable to incorporate NLO cross sections without any approximations in future analyses of 
parton densities. 

In the polarized case it is in general not at all clear whether a strategy based on correction 
factors Ki will work. Here, the parton densities are known with much less accuracy so far. It is 
therefore not possible to pre-calculate higher-order correction factors that one would be able to 
keep fixed throughout the fit, while using "fast" lowest-order expressions for the partonic cross 
sections. For instance, even though it is well known that for a sizable Ag the q + g — > 7 + q 
Compton subprocess is the dominant contributor to the transverse momentum distribution of 
a "prompt" photon in the kinematical region of interest, this is by no means the case if Ag 
happens to be small, in which case all other channels, even genuine NLO ones, may become 



equally important. In addition, the spin-dependent parton distributions, as well as the polarized 
partonic cross sections, may have zeros in the kinematical regions of interest, near which the 
predictions at lowest order and the next order will show marked differences. Therefore, even 
if the correction factors are updated at times during the fitting procedure, the convergence of 
the fit is not warranted. Conversely, if one updates the K t frequently, the fit will become too 
slow again. 

Clearly, in the polarized case, the goal must be to find a way of implementing efficiently, and 
without approximations, the exact NLO expression for any hadronic cross section such as the 
prompt photon cross section into the fitting procedure. As will be shown in the next Section, 
this can be achieved in a very simple and straightforward way by going to Mellin-n moment 
space. A technique of this sort was first used for the case of jet production in DIS as a means 
of extracting information about the unpolarized gluon density ||. The relevant generalization 
to hadron-hadron scattering, which is more involved and requires a "double Mellin transform" 
was recently provided in |J. However, || focuses on the formalism and the technical aspects 
of the Mellin transformation, rather than on its actual practicability, and the usefulness in a 
global QCD analysis has never been demonstrated. 

Before we demonstrate in some detail the potential of the Mellin technique in praxis for two 
examples relevant for future global analyses of polarized parton densities, which is the main 
thrust of this paper, we start off in the next Section by rederiving the required formalism in an 
easy and transparent way. In Section 3 we will consider first the semi-inclusive production of 
hadrons in polarized DIS as the simplest application of the Mellin technique. The n moments 
of the partonic cross sections can be taken analytically in this case. Due to the subsequent 
fragmentation of a final state parton into the observed hadron, semi-inclusive DIS (SIDIS) is 
sensitive to different flavor combinations than inclusive DIS data. It has also the advantage 



that we have already data at our disposal JR], [1JJ which can be analysed in terms of a possible 



flavor asymmetry Am — Ad of the light sea. As a second example we study the production 
of a prompt photon at high transverse momentum p? in pp collisions at RHIC in Section 4. 
Its sensitivity to the gluon distribution via the LO Compton subprocess, which, along with 
the cleanliness of the prompt photon signal, is the reason why this process will be the flagship 
measurement of Ag at RHIC 0. As a first case study for future global analyses we also carry 
out a "toy" analysis of DIS and projected prompt photon data to highlight the power of future 
RHIC pp data to pin down Ag. We briefly summarize the main results in Section 5. 



2 Hadronic Cross Sections and the Mellin Moment Tech- 
nique 

The factorization theorem [T2] ensures that in the presence of a hard scale in a reaction the cor- 
responding (spin-dependent) hadronic cross section can be written as a sum over "convolutions" 
of parton densities with partonic hard-scattering cross sections. The latter are perturbatively 
calculable and are specific to the reaction under consideration. The parton distributions, which 
for spin-dependent interactions contain the desired information on the nucleon's spin structure, 
depend on long-distance phenomena. However, they are universal: a single set of distributions 
for (anti-) quarks u,u,d,d,s,s, . . . and gluons g, predicts all data sets simultaneously. 

To be specific, for a general spin- dependent cross section in longitudinally polarized pp 
collisions, differential in a certain observable O and integrated over experimental bins in other 
kinematical variables T, one has 
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(1) 



= 5Z / dT dx a dx b dz c Af a (x a ,fi F )Af b (x b ,fi F )D^(z c ,fi' F ) 

} r*vr\ Kin / ™min / ~,min / _.min 

a,b,c t/ex P Dm Jx a ° x b J z c 

A A £ c 

x dQ jt (x a P A , x b P B , P H /z c , T, [i R , fi F , fi' F ) , 

where the arguments (++) and (H — ) in the first line of Eq. ([!]) refer to the helicities of the 
incoming hadrons A and B. The A/j are the spin-dependent parton distributions, defined as 

Afi(x, fi F ) = /+(x, up) - fr( x , (i F ) , (2) 

where ff (ff) denotes the number density of a parton-type fa with helicity '+' ('— ') in a 
proton with positive helicity, carrying the fraction x of the proton's momentum. The D^(z, fj,' F ) 
represent the unpolarized fragmentation functions. They parameterize the probability that a 
parton c fragments into the observed final state H, e.g., a charged pion, with momentum 
Ph — z Pc- For some observables, such as (di-)jets, there is no need for a fragmentation 
function in Eq. (jl]). 

The scales \i F and /j,' F are the factorization scales for initial and final state collinear singu- 
larities, respectively, and reflect the certain amount of arbitrariness in the separation of short- 
distance and long-distance physics embodied in Eq. ([!]). Even though the parton densities 
(fragmentation functions) cannot presently be derived from first principles, their dependence 
on fi F (fjL F ) is calculable perturbatively in terms of the "DGLAP" evolution equations [p~3| , al- 



lowing to relate their values at one scale to their values at any other fi F (fx' F ). The other scale, 



fiR, in Eq. (|l|) is the renormalization scale, introduced in the procedure of renormalizing the 
strong coupling constant. Finally, the sum in Eq. (p]) is over all contributing partonic channels 
a + b — > c + X, with dAa^ b the associated partonic cross section, defined in complete analogy 
with the first line of Eq. (|l|), the helicities now referring to partonic ones: 



dAa c ab = - 



dcr c ab (++) ~ dcr c J 



(3) 
As mentioned earlier, the dAa^ b are perturbative, that is, they have the expansion 

dMl„ = d^T + (£) dA»T + (^) 2 **#? » + . . . . (4) 

It should be noted that lepton-hadron reactions are also included in Eq. (|l|) by simply 
setting Afb(xb, up) — ^(1 — x b)- We will consider this example in some detail in Section 3 as 
it is the simplest application of the Mellin moment technique which we are going to advocate 
in the following as a straightforward tool to extract information about parton densities from a 
global QCD analysis. 

For the polarized parton distribution functions, the Mellin moments are defined as 

Af?(ji)= f dxx^Afifan). (5) 



It is well known [1 1] that the evolution equations for the the parton densities become particularly 
simple in Mellin-n space, since the convolutions occuring in the x-space equations factorize into 
simple products under moments. This allows for a straightforward analytic solution of the 
differential evolution equations, see, e.g., |L4|]. In fact, several of the NLO evolution codes used 
for parton density analyses in the unpolarized and polarized cases are set up in Mellin-n space. 
After evolving from one scale to another in moment space, the evolved parton distributions in 
Bjorken-rc space are recovered by an inverse Mellin transform, given by 



A/ i (x,/*) = ^ / dnx- n Af?(ri 

2m Jc n 



(6) 



where C n denotes a contour in the complex n plane that has an imaginary part ranging from 
— oo to oo and that intersects the real axis to the right of the rightmost poles of the A/"(/z). 
The evolution of the time-like fragmentation functions can be treated in a very similar way in 
Mellin space as well. 

The crucial, but simple, step in applying moment techniques to Eq. (|1|) is to express the 
Afi(x{, ftp) by their Mellin inverses in Eq. (Q) |J. One subsequently interchanges integrations 



and arrives at 
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b 



-n^ — m r)H I _, ,,/ 



dAa c 
x , Q ,y (xo.Pa, XbPs , Pfr/« c , r, ijlr, n F , n'p) (7) 

- J2 f dn f dmAf2(iM F )Afr(ii F )A&*( ni m,0,ii R ,ii F ). (8) 

One can now pre-calculate the quantities Aa^ b (n, m, O, [ar, fi F ), which do not depend at all on 
the parton distribution functions, prior to the fit for a specific set of the two Mellin variables 
n and m, for each contributing subprocess and in each experimental bin. Effectively, one has 
to compute the cross sections with complex "dummy" parton distribution functions x~ n x b m . 
We emphasize that all the tedious and time-consuming integrations are already dealt with 
in the calculation of the Acr^(n, m, 0,/i^, fi F ). We have included the integration over the 
fragmentation function D^ and the summation over the final state parton c in the definition of 
the pre-calculated quantities Aa^ b . This also implies that Aa^ b does not depend anymore on 
the choice for \J F apart from some residual dependence which is of higher order in a s . Usually 
the fragmentation functions are taken from "elsewhere", i.e., e + e~ data, rather than being 
fitted simultaneously with the parton densities. We note, however, that one can also replace 
D^ by their Mellin inverse according to Eq. (|J). In that case the pre-calculated quantities 
would depend on three Mellin variables. 

The double inverse Mellin transformation which finally links the parton distributions with 
the pre-calculated Aa^ b (n,m,0, /i#, n F ) of course still needs to be performed in each step of 
the fitting procedure. However, the integrations over n and m in Eq. @ are extremely fast to 
perform by choosing the values for n, m in Aa^(n, m, O, /i#, fi F ) on the contours C n , C m simply 
as the supports for a Gaussian integration. The point here is that the integrand in n and m falls 
off very rapidly as \n\ and \m\ increase along the contour, for two reasons: first, each parton 
distribution function is expected to fall off at least as a power (1 — a;) 3 at large x, which in 
moment space converts into a fall-off of ~ l/n A or higher. Secondly, we may choose contours in 
moment space that are bent by an angle a — n/2 with respect to the vertical direction; a possible 
choice is shown in Fig. |l|. Then, for large \n\ and \m\, n and m will acquire large negative real 
parts, so that (x a )~ n and (x&) _m decrease exponentially along the respective contours. This 
helps for the numerical convergence of the calculation of the Aa^ b (n,m,0,p,R,p, F ) an d also 
gives them a rapid fall-off at large arguments. We note that no new poles in n and m, beyond 
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those already present in the moments of the parton distribution functions, are introduced by 
the Acr^(n, m, O, p, R , p, F ) 0. 

We note that if one wishes to integrate also over an experimental bin in O in Eq. ((?]), a 
potential complication arises if the hard scale \ip in the parton distribution functions depends 
explicitly on O. This makes it impossible to straightforwardly include the O integration in 
the pre-calculation of the Aa^ b . A typical example for O, which often appears in practice, 
is the transverse momentum p? of an observed jet, hadron, or prompt photon. In this case, 
the O dependence of p,p is, however, not a serious limitation ||: the logarithmic dependence 
of the parton densities on \xf is much weaker than the overall pr dependence of the cross 
section. Therefore, it is always possible to choose a bin-average of px as the scale in the parton 
densities. Alternatively, one could choose not to include the pt integration in the Aa^ b and to 
construct grids of somewhat larger size, taken at a small number of support points for a simple 
Gaussian integration over the pt bin. A further possibility |§ is to absorb also the evolution 
of the parton densities from their initial scale po to \xf into the Aa^ b , which in moment space 



simply enters in the form of exponentials involving the anomalous dimensions, see, e.g., [[14 
This procedure, which is somewhat more involved, would eliminate any complication related 
to up ~ 0(px)- Anyway, the experiments will usually quote results for the pr- differential cross 
section at the pr-average over the bin, which of course is exactly what we have considered in 
Eqs. (0), (H). In the latter case it is also easily possible to organize the grids in such a way that 
the renormalization and/or factorization scales can be varied during the fit, by simply taking 
the (logarithmic) dependence on p R>F and the strong coupling a s (p R ) out of the partonic cross 
sections beforehand. 

As a technical sidestep, we give an explicit expression for the double inverse transform in 
Eq. (H) for the contours depicted in Fig. |l|. To this end, we parameterize the various segments 
in Fig. | by 

n = c n + u n e ±ia and m = c m + u m e ±ia , (9) 

where u ntm G [0, oo] and the sign of a has to be chosen appropriately for the branches of the 
contours. We then find 

-VRe / du n du m Af:ifi F ) (10) 

71 „ i, l Jo Jo 



dO 2.. 

a,b 



x|e aa A/rWA^(n,m,0, to ^) - (A^))* Acr*(n,m*, 0^ R ^ F ) 

where the asterisk denotes the complex conjugate, and where we have made use of A/™* (/xp) = 
(A/ b m (/ip))*, since the Afi(x, p) are real functions. This identity also implies that there is no 




Figure 1: Contours in complex Mellin-n,m spaces for the calculation of the double Mellin inverse 
in Eq. (|). r n and r m denote the rightmost poles of the integrand in n and m, respectively, and the 
c, the intersections with the real axis. 



need to separately compute the moments of the parton densities at the complex conjugate values 
n*,m*, which has a further positive effect on the computing time required for performing the 
Mellin inverses. In addition, there is no need to provide separate grids for Aa^ b (n*, m, O, //r, fip) 
and Aa^ b {n*,m*,0,fi R ,fi F ). 

Before proceeding, we reemphasize that the idea outlined above of reverting to Mellin mo- 
ment space in the implementation of any higher-order cross section into parton density fits is 
not entirely new, but was first developed in Refs. [R 0]. The example considered in || was jet 
production in DIS, which offers the simplification of being only linear in the parton distribution 
functions. There is a difference between our approach and that of Ref. || in practical terms: 
in the language of our example in Eq. (|7j), Ref. || would insert a factor (x™ m ) Tl (x™ m ) m in 
the integrands for the x a and Xb integrations, while undoing this operation through a factor 
(x™ m ) _n (x™ m ) _m in the n and m integrands. Even though obviously equivalent mathemati- 
cally, the disadvantage of this procedure is that the resulting factors (x™ m /x a ) n , (x™ m / Xb) m in 



the x a , %b integrands will now grow exponentially along the contours in Mellin space, making 
it numerically much more cumbersome || to perform the x a and x^ integrations yielding the 
AcH^(n, m, 0,hr, fi F )- As a matter of fact, the actual extension of the method of Ref. || to 
the case of hadronic collisions involving bilinear combinations of parton distributions appears 
difficult. The generalization of || to hadron-hadron scattering, without the shortcomings men- 
tioned above, was first provided in ||. The difference between our organization of the expression 
in Eq. (||) and Ref. || is the choice of the contour. Reference fully exploits the freedom in 
deforming the contours for the inverse Mellin transform and constructs a "surface of steepest 
descent" which in principle has the best numerical convergence properties but is difficult to 
parameterize. Instead we stick to the simple contours in Fig. p] which, as we will show below, 
turn out to be sufficient to obtain numerical agreement between Eqs. (|1|) and (P) of far better 
than 1% for all applications we are going to consider. We should also note that in |J the 
usefulness of the Mellin transform method was not demonstrated in practice. 

3 Semi-Inclusive Deep-Inelastic Scattering 



As a first application for the Mellin transform technique outlined in the previous Section we 
consider the semi-inclusive production of a hadron H in DIS. SIDIS starts at the Born level with 
the LO reaction 7*g — > q. The NLO 0(a s ) corrections also comprise the processes 7*g — > qg 
and 7*g — > qq and have been calculated in the spin-dependent case in the MS scheme in [|15[| . 
In each case one of the final state partons subsequently fragments into the observed hadron H. 
As in the fully inclusive case, the expression for the cross section is given by a single structure 
function gi(x, z, Q): 



dAa 



H 



Aira 2 



dx dy dz Q 2 



{2-y)g?(x,z,Q). 



To NLO in a s , gf can be written as [p~5] , [16] 
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Ag(|,^)ACW(a 
A 9 (l,, F )AC S (^/f)Df(y F 



(11) 



(12) 



with x and y denoting the usual DIS scaling variables (Q 2 = —q 2 = xyS), and where [|17], [18 



z = pn ■ Pn/Pn ■ Q- Eq. (|12"D and the variable z only apply to hadron production in the current 
fragmentation region characterized by positive values for the Feynman variable xf- All NLO 



MS partonic coefficient functions AC}- are collected in App. C of |T3|. They are non-trivial 



functions of x and z such that the x and z dependences of the cross section do not factorize into 
separate functions. Therefore the inclusion of the NLO corrections seems to be indispensable 
for a reliable extraction of parton densities from SIDIS. 



Due to the double convolutions appearing in Eq. ([E|) and the fact that the coefficient func- 
tions contain mathematical distributions as x — > 1 and/or z — > 1, the direct use of Eq. (|12|) 
in a global analysis of parton densities is rather time consuming and awkward, though not 
impossible |B5| since the partonic coefficient functions are still fairly simple. SIDIS has the ad- 
vantage, however, that the Mellin moments in x and z can be taken completely analytically for 
the partonic coefficient functions in Eq. (|T2|) . In doing so, the double convolutions in Eq. (|L2]) 
reduce to simple multiplications ~ A/ n (/i^)AC|- m (^, *-q)DV % ([i f ) and all distributions be- 



come ordinary functions of the moment variables. The ACJ- , defined by 






dxx n ~ l / dzz m ~ x 
Jo 



AC v (X ' Z ' Q> Q } ' 



(13) 



W/ 



are straightforwardly determined from the expressions for the corresponding AC} ,- (x, z, ^-, -§■) 
in App. C of [[TJ|] and read: 
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(15) 
(16) 
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where C F = 4/3, Tr = 1/2, and 



n 



(17) 



0=1 



For completeness we give also the Mellin moments for the corresponding unpolarized coeffi- 
cient functions C[ i J 71171 and C L i 'J im relevant for the structure functions Ff 1 and F^ , respectively. 
Using again Eq. (|i~3"D and the x, z space expressions in App. C in |L5| one finds 
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(21) 
(22) 
(23) 



'(n+l)(n + 2) 

We note that the usefulness of taking double Mellin moments for unpolarized SIDIS was first 
pointed out, though not further pursued, in [I7|. In the polarized case Mellin- n moments of 
the semi-inclusive cross section at fixed z have been recently considered in 



Having available the coefficient functions in Mellin moment space one can evaluate the 
desired SIDIS structure function g^ in a fast way by a double inverse Mellin transform as 
discussed in Section 3. One further ingredient required is the evolution of the moments of 
the fragmentation functions DV l (fJ /F ) which proceeds along very similar lines as for the parton 



densities. Below we will use the recent NLO analysis of [21] which can be applied down to the 



Q values required for the available spin-dependent SIDIS fixed-target data [JT0|, [Tljl . It should 
be noted that the Mellin approach allows in principle a simultaneous fit of parton densities and 
fragmentation functions in SIDIS at no extra "costs" . 

The experimentally relevant quantity is the so-called spin asymmetry, defined as the ra- 
tio of the polarized and unpolarized SIDIS structure functions, g^(x,z,Q) and F^(x,z,Q), 
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Figure 2: Comparison of a fit to DIS and SIDIS data in NLO QCD (see text) with the measured 
SIDIS spin asymmetries A x ' for the production of positively or negatively charged hadrons (±) off 



different targets N |L0], |Tl| 



respectively, 



A h (t 7 (1) - 9l ^ifl^j 



(24) 



Due to limited statistics all presently available results for Eq. (|24] ) are integrated over the 
entire z range accessible experimentally (z > 0.2) |I| 11[. To facilitate the comparison with 
these data it is more convenient to define an "effective" coefficient function AC rather than 
using the double moments and integrating afterwards over z. The AC already incorporate the 
z integration and can be easily pre-calculated once prior to the fit. They are defined by 



ACf^) ee f dz [ dmz- m AC^ nm (^^)DT^' F 
V Jz min Jc m V V 



(25) 



(in LO one has AC/- = 1) and can be used in a similar way as the usual fully inclusive DIS 
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Figure 3: Resulting flavor asymmetry of the light sea, Aw — Ad, at the input scale /i from our 
combined fit to DIS and SIDIS data. Also shown are model predictions taken from @ (dotted line) 
and M, 123 (dashed line). 



coefficient functions. This makes the numerical evaluation extremely fast: 100 calculations of 
the SIDIS cross section in NLO take only about 1 second on a standard workstation. Clearly, 
SIDIS data can be as easily incorporated in a global QCD analysis as DIS data. 

In Fig. we compare the result of a NLO fit to all available data for DIS and SIDIS spin 
asymmetries with data for A l ' for positively or negatively charged hadrons H ± and different 



targets iV [pll ITTJ. Regarding the details of the analysis, we stay in the framework of the 



"standard" fit of [22|, but allow for an SU(2) breaking of the light sea by introducing a function 
/su(2) j 



Au' (x, no) = Au(x,Hq) - /su(2)(z,A*o) 
Ad'(x, no) = Ad(x, fi ) + / SU (2) {x, /i ) 



Au'{x, Ho) = Au(x, Ho) + /su(2) (ar, Ho) , 
Act (x, Ho) = Ad(x,Ho) - fsv(2)(x,Ho) , (26) 



such that all quark combinations measured in inclusive DIS remain unchanged, Aq' + Aq' = 
Aq + Ag, but Am' — Ad 1 = 2/su(2) • We choose a "minimal" ansatz for /su(2) with 3 additional 
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parameters 

fsu(2)(x,ti ) = Nx a (l-x)P (27) 



where /io — 0.6 GeV is the initial scale for the evolution in p2| . We choose the renormalization 
and factorization scales //r = Hf = h'f = Q- 

The resulting asymmetry of the light sea at the input scale /^o is shown in Fig. 0. For 
comparison we also show model predictions for Am — Ad from || and 0, p2[ . It turns out that 
the flavor asymmetry obtained in our analysis is much less pronounced than predicted in most 
models. It has to be stressed, however, that the change in the total \ 2 f° r a U SIDIS data is 
less than one unit if one chooses an SU(3) symmetric sea, the model calculations [^ ^, |2!| or 
our fit result. Thus one has to conclude that present SIDIS data are not precise enough to 
distinguish between different results for Am — Ad and that one has to wait for new SIDIS data 
from HERMES and, in particular, for results on W ± boson production at RHIC ||. Similar 



conclusions have been reached in the analysis of |T9| 



4 Prompt Photon Production at RHIC 

To give an example for the Mellin technique in hadron-hadron collisions, we study the produc- 
tion of a prompt photon in pp collisions at RHIC. In this case dAa^l in Eq. (|) starts at LO 
with the reactions q + q — > 7 + g and q + g — > 7 + q, the latter channel being sensitive to the 
polarized gluon distribution. The NLO corrections, dAa^l , are also available [23]. The NLO 



x-space expressions are rather lengthy and complicated, and Mellin moments cannot be taken 
analytically anymore. Nevertheless, as we shall see below, it is in the analysis of hadron-hadron 
collision data where the Mellin moment technique exhibits its full potential and usefulness. 

To be specific, the transverse momentum (pt) distribution of a prompt photon in pp colli- 
sions at a center-of-mass energy yS, integrated over a certain experimental bin in pseudora- 
pidity r] [i.e., "O = pr" and "T = 77" in Eq. (JlJ) ] , is given by 



—, = y^ / dr] / dXa I dx b Af a (x a , fi F ) Af b (x b , fi F ) 

UPT ab Jri-bin Jx™ in Jx™ in 

X- — j^(x a P A , x b P B ,p T ,r], n R ,fx F ) , (28) 

dprdri 

where x™ in = XTe v /(2 — i r e"' ) ) and x™ in = x a XTe~ v /(2x a — xtQ v ) with x F = 2p T /\/S. For 
our case study, we analyse the polarized prompt photon cross section in Eq. (|28| ) at NLO. 
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Figure 4: Comparison of the results based on the Mellin-technique in Eqs. (0) and (|8]) to those of 
Eq. (|^) for various sizes of the grid in n,m. 



The associated spin-asymmetry, denned as the ratio of the polarized and the unpolarized cross 
sections, 

(29) 



A 1 



dAa 1 /dp? 



da"'/dpT 

will soon be measured at RHIC in collisions of longitudinally polarized protons and, as men- 
tioned above, will be a key process for measuring Ag. We use v^S* = 200 GeV and look at 
the cross section as a function of the photon's transverse momentum pt for five values of px 
which will be experimentally accessible at RHIC, px = [12.5, 17.5, 22.5, 27.5, 32.5] GeV. We 
average over \q\ < 0.35 in pseudorapidity. As in experiment ||, we impose an isolation cut 
on the photon, for which we choose the isolation proposed in [Q with parameters R = 0.4, 
e = 1. A positive feature of this isolation criterion is the absence of a fragmentation contribu- 
tion to prompt photon production, hence we can drop the z c integration and the fragmentation 
function D^ in Eq. ([I]). We choose the renormalization and factorization scales fiR = Hf = Pt- 

Our first goal here is to show that the method based on Eqs. (0) and (H) actually works also 
for the more complicated case of hadron-hadron collisions and correctly reproduces the result 
obtained within the direct, but "slow", calculation via Eq. (f28|). Also, we need to establish an 
optimal size of the grids that yields excellent accuracy but is still calculable in, say, a few hours 
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of CPU time on a standard workstation. Fig. |] compares the results based on Eqs. (^) and 
(§), referred to as "Mellin technique", to those of Eq. fl28|) , for various sizes of the grid in n, m. 
Here we have used again the polarized parton densities of 2!| ("standard" set). For a more 



detailed comparison, we split up the contributions to the NLO prompt photon cross section 
into three parts, associated with the reactions q + q — > 7 + X and q + g —>■ 7 + X that are 
already present at Born level, and all other processes that arise only at NLO. One notices that 
in each case already a grid size of 64 x 64 values yields excellent accuracy. Even a 56 x 56 grid 
is acceptable apart from a minor deviation occuring for qq scattering in the vicinty of a zero 
in the partonic cross section. We have checked that the results in Fig. |] do not depend on the 
actual choice of parton densities. 

The crucial asset of the Mellin method is the speed at which one can calculate the full 
hadronic cross section, once the grids AaJ b (n,m,p T , fi R , fi F ) have been pre-calculated. For the 
64 x 64 grid, we found that 1000 evaluations of the full NLO prompt photon cross section 
take only about 10-15 seconds on a standard workstation. Note that this number includes the 
evolution (in moment space) of the parton distributions from their input scale to the scale pr 
relevant to this case. Clearly, an implementation into a full parton density fitting procedure is 
now readily possible. 

To give an example, we finally perform a "toy" global analysis of the available data on 
polarized DIS [l|] and of fictitious data on prompt photon production at RHIC |J, which we 
project by simply calculating A i LL in Eq. (|29"D to NLO using the sets of polarized and unpolarized 
parton distributions of [^2j and 0], respectively. For an estimate of the anticipated la errors 
on the "data" for A1 L , we use the numbers reported in ||. We subsequently apply a random 
Gaussian shift of the pseudo-data, allowing them to vary within la. The "data", as well as the 



underlying theoretical calculation of A\ L based on the spin- dependent parton densities of p2 | 
(solid line), are shown in the left panel of Fig. |[ 

Next, we perform a large number of fits to the full, DIS plus projected prompt photon, data 
set. We simultaneously fit all polarized parton densities, (anti)quarks and gluons, choosing the 
distributions of |22[] as the input for the Ag, Ag, but using randomly chosen values for the 
parameters in the ansatz for the polarized gluon distribution at the input scale /i . Regarding 
the details of the evolution, we stay again within the setup of \T1\ , but we choose a more flexible 
ansatz for the polarized gluon density, 

Ag(x,fio) = Nx a (l-x) f3 (l + 1 x)g(x,fio), (30) 

which also allows for a zero in the x-shape of Ao. o(x,/i ) is the unpolarized gluon density [|7| 



at the input scale of [22| . Note that the functional form for the polarized gluon density of 
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Figure 5: (a): Generated pseudo data for A\ L based on a calculation using the spin-dependent 
parton densities of [j22[] (solid line). The shaded band corresponds to the results of a large number 
of combined fits to DIS and A\ L data (see text), (b): Variations (shaded area) of the polarized 
gluon densities obtained in the combined fits, with respect to Ag of |^5], for fx F = 10 GeV (see 
text). Also shown are two extreme gluon densities (dotted lines) which give an excellent description 
of the polarized DIS data only. 



used for generating our pseudo-data, is included in Eq. ( |30| ) for N = 1.419, a 
0.15, 7 = 0. Each fit takes only about 10 -r- 20 minutes. 



1.43, (3 



Ideally, thanks to the strong sensitivity of the prompt photon reaction to Ag, the gluon 
density in each fit should return close to the function we assumed when calculating the fictitious 
prompt photon "data" , in the region of x probed by the data. Indeed, as shown in Fig. ||(b), this 
happens. The shaded band illustrates the deviations of the gluon densities obtained from the 
global fits to the "reference Ag" p2[ used in generating the pseudo-data. It should be stressed 
that only those fits are admitted to the band that give a good simultaneous description of the 
DIS and A 1 LL data. Here we have tolerated a maximum increase of the total \ 2 by up to four 
units from its minimum value. The shaded area in Fig. [5](a) shows the corresponding variations 
in A\ L . 

As is expected, all gluon densities are rather tightly constrained in the x-region dominantly 
probed by the prompt photon data. This is true in particular at x ~ 0.15, as a result of the 
most precise data point for A\ L at pt = 12.5 GeV. We note that one can also easily include 
the SIDIS data discussed in Section 3 into the global analysis without any significant increase 
of computing time for each fit. However, so far these data have no impact on our results. To 
illustrate our present ignorance of Ag, Fig. |5|(b) shows also two extreme gluon densities with 
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first moments Ag 1 (fio) = ±0.8 (dotted lines), which are both in perfect agreement with all 
presently available DIS data. The corresponding predictions for A y LL for these two sets are given 
in Fig. H](a). It should be noted that future measurements of K* LL at RHIC at \[S = 500 GeV 
and for similar px values of the prompt photon would further reduce the uncertainties on Ag 
in the x-region between 0.05 and 0.1. Although our analysis still contains a certain bias by 
choosing only the framework of JEJ for the fits as well as by our choice of what \ 2 values are 
still tolerable, it clearly outlines the potential and importance of upcoming measurements of 
A\ L at RHIC for improving our understanding of the spin structure of the nucleon, in particular 
of its spin-dependent gluon density. 



5 Conclusions 

To conclude, we have presented and applied a powerful technique for implementing in a fast 
way, and without any approximations, higher-order calculations of partonic cross sections into 
global analyses of parton distribution functions. We have demonstrated that the approach 
works in practice for two examples: SIDIS and prompt photon production in pp collisions. 
In the first case it was possible to perform the Mellin transform analytically and we have 
provided all necessary technical details for future analyses of polarized and unpolarized SIDIS 
data. For polarized prompt photon production we have presented a case study for a future 
global analysis based on fictitious data. The Mellin transform method is certainly applicable to 
any other reaction of interest, and it could equally well be an improvement also in any global 
analysis of unpolarized parton distributions. 
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